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Abstract 

We present a calculation of the next-to- leading order QCD corrections, with one-scale phase space 
slicing method, to single top quark production and decay process pp, pp ^ th + X ^ hivh + X aX 



hadron colliders. Using the helicity amplitude method, the angular correlation of the final state 
D , partons and the spin correlation of the top quark are preserved. The effect of the top quark width 

is also examined. 
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I. INTRODUCTION 



At hadron colliders, the top quarks (t) are predominantly produced in pairs via the strong 
interaction process qq, gg tt. Though it is possible to study the decay branching ratios 
of the top quark in tt pairs, to test the coupling of top quark with bottom quark (b) and 
W gauge boson in hadron collisions, it is best to study the single-top quark production. 
Compared to the top quark pair production, produced by the interaction of the Quantum 
Chromodynamics (QCD), the single top quark productions are through the electroweak 
interaction connecting top quark to the down-type quarks, with amplitudes proportional 
to the Cabibbo-Kabayaashi-Maskawa (CKM) matrix elements. Due to the nature of left- 
handed charged weak current interaction, the top quark produced via single-top processes 
is highly polarized. Furthermore, top quark will decay via weak interaction before it has 
a chance to form a hadron, so its polarization property can be studied from the angular 
distributions of its decay particles. Hence, measuring the production rate of the single-top 
event can directly probe the electroweak properties of the top quark. For example, it can be 
used to measure the CKM matrix element Vtb and to test the V—A structure of the top quark 



charged- current weak interaction, or to probe CP violation effects j^Q,^]- Besides of playing 
the role as a test of the Standard Model (SM), the precision measurement of the single top 
quark events has additional importance in searching for new physics, because the charged- 
current top quark coupling (W-t-b) might be particularly sensitive to certain new physics 
via new weak interactions or via loop effects, and new production mechanism might also 



contribute to the single top event rate 



aertects, ana new proauctic 



mi gnt also 



Furthermore, the single-top event is also an important background to the search of Higgs 



boson 
search 



qq 



WH with H bb) at the Tevatron 



0, H Q 



and other new physics 



24|. 



Because of the unique features of the single top quark physics, it has been extensively 
studied in the literature 

Q, HQ, HQ, HH0HHHHHHH 

Q,HQQQ,Q. There are three separate single top quark production processes of 
interest at the hadron collider, which may be characterized by the virtuality of the W boson 
(with four momentum q) in the processes. The s-channel process qq' — ^ W* tb via a 
virtual s-channel W boson involves a timelike W boson, > (m^ + m^)^, the t-channel 
process qb q't (including q'b qt, also referred as ly-gluon fusion) involves a spacelike 



9. 



W boson, < 0, and the tW associated production process bg — tW"" involves an on-shell 
W boson, = m^. Therefore, these three single top quark production mechanisms probe 
the charged-current interaction in different regions and are thus complementary to each 
other. Furthermore, they are sensitive to different new physics effects 4^, and should be 
separately studied. 

To improve the theoretical prediction on the single-top production rate, a next-to-leading- 
order (NLO) correction, at the order of a^, for the s- and t-channel processes has been carried 



out in Refs 



decay 



32 



39|. This is similar to the study of the 0{as) correction to the top quark 



. Although the above studies provide the inclusive rate for single-top production. 



they cannot predict the event topology of the single-top event, which is crucial to confront 
the theory with experimental data in w. 



such an event. For that, Refs. 



47 



/■hich some kinematical cuts are necessary to detect 
48| have calculated the differential cross section for 
on-shell single top quark production. However, NLO corrections to the top quark decay 
process were not included, nor the effects of the top quark width were considered. Since 
the top quark production and decay do not occur in isolation from each other, a theoretical 
study that includes both kinds of corrections is needed. A complete NLO calculation should 
include contribution from the production and the decay of the top quark, and the angular 
correlation among the final state particles should be calculated to analyze the polarization 
of the top quark. The 0{as) corrections to kinematic distributions may depend on the 
kinematic cuts and on the jet algorithm that must be implemented in the experiments. 
Therefore, it is necessary to obtain a fully differential calculation that can be used to study 
the kinematics of the final state particles. 

In this theory paper, we present a NLO QCD calculation with the one-scale phase space 
slicing method, which treats consistently 0{as) corrections to both the production and the 
decay of the top quark in single top events. Our approach can give not only the inclusive 
total cross section, but also the various kinematical distributions of the final state particles, 
and provide a study on the top quark polarization at the NLO. Furthermore, since realistic 
kinematical cuts can be applied, our approach allows the experimentalists to compare their 
results directly with the theoretical predictions. In our study, we assume in all cases leptonic 
decays of the W boson (for the sake of definiteness, we shall consider — > e^v] the lepton 
mass effects will be neglected throughout this paper). The phenomenological discussions 
will be given in our sequential papers 



The rest of this paper is organized as follows. In Sec. II, we outline the method of our 
calculation. In Sec. Ill, we present the Born level helicity amplitudes of the production 
and decay of single top quark. In Sec IV, we present the NLO helicity amplitudes of the 
production and decay of single top quark. The effective form factor approach is adopted 
in the calculation to generalize the application of our formalism to, for example, studying 
new physics effects. In Sec. V, we use the phase space slicing (PSS) method to calculate 
the effective form factors. To regularize divergencies in the calculation that involves the 
75 matrix, both the dimensional regularization (DREG) [s^ and the dimensional reduction 
(DRED) 1^ schemes are examined and their difference is shown in each individual form 
factor. In Sec. VI, we show how to assemble all the components discussed above to enumerate 
the NLO differential cross section of the single top quark. Finally, we give our conclusions 
in Sec. VII. 



II. OUTLINE OF THE CALCULATION 



In this section we outline the method of our calculation whose details shall be presented 
in the following sections. 



A. Narrow width approximation 



In this work, the narrow width approximation (NWA) is used to study the production 
and decay of single top quark, in which the 0{as) corrections can be unambiguously assigned 
to either the single top quark production process or the top quark decay process. A finite 
top width will result in a new type of virtual NLO Feynman diagram in which a gluon line 
is connected from the anti-bottom quark (of the production process) to the bottom quark 
(of the decay process). Moreover, there will also be interference between the gluons emitted 
in the production and the gluon emitted in the decay if the effects of finite top quark width 
is considered. Those effects are nonfactorizable, which are similar to the effects of QED 
radiative corrections to the scattering process e~^e~ — > W^W~ — > 4/. It was shown that the 



nonfactorizable effects are small as long as the process is not near the threshold 



This provides the motivation of using the NWA for this kind of calculation 



3, 



53, Q- 



The single top quark can be produced through s-channel and t-channel processes, as 



A 



shown in Fig. ^a) and Fig. ^b), respectively. Using the NWA, we decompose the Born 
level processes, depicted in Fig. ^as indicated by symbol (E), into two parts: the top quark 
production and its sequential decay, where both the production and decay matrices are 
separately gauge invariant. Making use of the polarization information of the top quark, we 
can apply the NWA to correlate the top quark production with the top quark decay processes 
by replacing the numerator of the top quark propagator {^t + f^t) by '^x^=±u''^^{t)u^^{t). 
Here u^^{t) is the Dirac spinor of the top quark with helicity Aj, where \t = + oi — for a 
right-handed or left-handed top quark, respectively. Therefore, the scattering amplitude of 
the single top quark production and decay processes can be written as [J] 

M=Y^ M'^^%\t)MP''°\\t), (1) 
\t=± 

where M.{\t) is the helicity amplitude and is the helicity eigenvalue of the single top 
quark produced in the intermediate state. The matrix element squared can be written as 
the product of the production part and the decay part in the density matrix formalism: 

\M\'= A.AiS,„,j, (2) 

\t,K=± 

where 

^XuK = MlJX,)MUK), (3) 
^x.,K = MlJX,)M,UK)- (4) 

In addition to the matrix elements, the phase space of the single top quark processes can 
also be factorized into the top quark production and decay for an on-shell top quark in NWA 
by writing the denominator of the top quark propagator as 



/ 



When the matrix element is calculated using the fixed rrit value, it is the usual NWA method. 
In this case, the invariant mass of the top quark decay particles will be equal to rrit (a fixed 
value) for all events. Reconstructing the top quark invariant mass from its decay particles 
is an important experimental task at the Tevatron and the LHC, it is desirable to have a 
theory calculation that would produce the invariant mass distribution of the reconstructed 
top quark mass with a Breit-Wigner resonance shape to reflect the non-vanishing decay 




[a) (6) 

Figure 1: Feynman diagrams of the Born level contribution to the production and decay of single 
top quark, (a) s-channel (b) t-channel 

width of the top quark (for being an unstable resonance). For that, we introduce the 
"modified NWA" method in our numerical calculation in which we generate a Breit-Wigner 
distribution for the top quark invariant mass in the phase space generator and then calculate 
the squared matrix element according to Eq. (0) with rrit being the invariant mass generated 
by the phase space generator on the event-by-event basis. In the limit that the total decay 
width of the top quark approaches to zero (i.e., much smaller than the top quark mass), the 
production and the decay of top quark can be factorized. Therefore, the S-matrix element 
for the production and the decay processes are separately gauge invariant with any value of 
top quark invariant mass. We find that the total event rate and the distributions of various 
kinematics variables (except the distribution of the reconstructed top quark invariant mass) 
calculated using the "modified NWA" method agree well with that using the NWA method. 
In the NWA method, the reconstructed top quark invariant mass distribution is a delta- 
function, i.e., taking a fixed value, while in the "modified NWA" method, it is almost a 
Breit-Wigner distribution. The reason that the "modified NWA" method does not generate 
a perfect Breit-Wigner shape in the distribution of the top quark invariant mass is because 
the initial state parton luminosities (predominantly due to valence quarks) for the s-channel 
single-top process drop rapidly at the relevant Bjorken-x range, where (x) ~ ^ ~ 0.1. 

B. Phase space slicing method 

When calculating NLO QCD corrections, one generally encounters both ultraviolet (UV) 
and infrared (IR) (soft and collinear) divergencies. The former divergencies can be removed 
by proper renormalization of couplings and wave functions. We don't need to renormalize 



the couplings in our calculation because the Born level couplings do not involve QCD in- 
teraction. In order to handle the latter divergencies, one has to consider both virtual and 
real corrections. The soft divergencies will cancel according to the Kinoshita-Lee-Nauenberg 



(KLN) theorem [58|,|59l, but some collinear divergencies will remain uncanceled. In the case 
of considering the initial state partons, one needs to absorb additional collinear divergencies 
to define the NLO parton distribution function (PDF) of the initial state partons. After that, 
all the infrared-safe observables will be free of any singularities. To calculate the inclusive 
production rate, one can use the dimensional regularization scheme to regularize divergen- 
cies and adopt the modified minimal subtraction (MS) factorization scheme to obtain the 
total rate. However, owing to the complicated phase space for multi-parton configurations, 
analytic calculations are in practice impossible for all but the simplest quantities. During 
the last few years, effective numerical computational techniques have been developed to cal- 
culate the fully differential cross section to NLO and above. There are, broadly speaking, 
two types of algorithm used for NLO calculations: the phase space slicing method, and the 
dipole subtraction method 

0, 0, 0, 0, Q, 0, Q, 0, S Q- 

In this study, we use 

the phase space slicing method (PSS) with one cutoff scale for which the universal crossing 
functions have been derived in Ref. j63|. The advantage of this method is that, after calcu- 
lating the effective matrix elements with all the partons in the final state, we can use the 
generalized crossing property of the NLO matrix elements to calculate the corresponding 
s-channel or t-channel matrix elements numerically without requiring any further effort. The 
validity of this method is due to the property that both the phase space and matrix element 
of the initial and final state collinear radiation processes can be simultaneously factorized. 
Below, we briefly review the general formalism of the NLO calculation in PSS method with 
one cutoff scale. 

The phase space slicing method with one cutoff scale introduces an unphysical parameter 
Smin to Separate the real emission correction phase space into two regions: (1) the resolved 
region in which the amplitude has no divergencies and can be integrated numerically by 
Monte Carlo method; (2) the unresolved region in which the amplitude contains all the soft 
and collinear divergencies and can be integrated out analytically. It should be emphasized 
that the notion of resolved/unresolved partons is unrelated to the physical jet resolution 
criterium or to any other relevant physical scale. In the massless case, a convenient definition 
of the resolved region is given by the requirement Sij > Smin for all invariants Sij = (pi+Pj)'^, 



7 



where pi and pj are the 4-momenta of partons i and j, respectively. For the massive quarks, 
we follow the definition in Ref. [3| to account for masses, but still use the terminology 
"resolved" and "unresolved" partons. In the regions with unresolved partons, soft and 
collinear approximations of the matrix elements, which hold exactly in the limit s^m 0, 
are used. The necessary integrations over the soft and collinear regions of phase space 
can then be carried out analytically in d = 4 — 2e space-time dimensions. One can thus 
isolate all the poles in e and perform the cancellation of the IR singularities between the 
real and virtual contributions and absorb the leftover singularities into the parton structure 
functions via the factorization procedure. After the above procedure, one takes the limit 
e ^ 0. The contribution from the sum of virtual and unresolved region corrections is finite 
but Smin dependent. Since the parameter s^m is introduced in the theoretical calculation for 
technical reasons only and is unrelated to any physical quantity, the sum of all contributions 
(virtual, unresolved and resolved corrections) must not depend on Smin- We note that the 
phase space slicing method is only valid in the limit that Smin is small enough so that a given 
jet finding algorithm (or any infrared-safe observable) can be consistently defined even after 
including the experimental cuts. 

In general, the conventional calculation of the NLO differential cross section for a process 
with initial state hadrons Hi and H2 can be written as 

d^mS, = J2j ^^1^^2/fHa;i,/XF)/f^(x2,/iF)t^^5^(xiX2,/iij), (6) 
a,b 

where a, b denote parton flavors and Xi, X2 are parton momentum fractions, f^^x^^p) 
is the usual NLO parton distribution function with the mass factorization scale and 
da^f^^*-^ {xi, X2, fin) is the NLO hard scattering differential cross section with the renormal- 
ization scale fiR. The diagrammatic demonstration of Eq. (jH)) is shown in the upper part of 
Fig. 121 

Unlike the conventional calculation method, the PSS method with one cutoff scale will 
first cross the initial state partons into the final state, including the virtual corrections and 
unresolved real emission corrections. For example, to calculate the s-channel single top 
quark production at the NLO, we first calculate the radiative corrections to W* QQ'ig), 
as shown in the lower part of the Fig.|21 in which we split the phase space of the real emission 
corrections into the unresolved and resolved region. After we integrate out the unresolved 
phase space region, the net contribution of the virtual corrections and the real emission 

8 
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Figure 2: Illustration of the PSS method with one cutoff scale to describe the processes with initial 
state massless quarks. Here, only half of the real emission diagrams is shown. In this paper, we 
assign the particle's momentum such that the initial state particle's momentum is incoming to the 
vertex while the final state particle's momentum is outgoing. 



corrections in the unresolved phase space is finite but theoretical cutoff (smin) dependent, 
and it can be written as a form factor (denoted by the box in Fig. ^ of the Born level vertex. 

Secondly, we take the already calculated effective matrix elements with all the partons 
in the final state and use the universal "crossing function", which is the generalization of 
the crossing property of the LO matrix elements to NLO, to calculate the corresponding 
matrix elements numerically. Once we cross the needed partons back to the initial state, 
the contributions from the unresolved collinear phase space regions are different from those 
with all the partons in the final state. These differences are included into the definition of 
the crossing function (including the mass factorization effects), as shown in the middle part 
of Fig. El In this paper, we only present the explicit expressions of the crossing function. 



For t 
Ref. 



le definition and detailed derivation of the crossing function, we refer the reader to 
6^. After applying the mass factorization in a chosen scheme, the Smin dependent 



crossing functions for an initial state parton a, which participates in the hard scattering 



processes, can be written in the form: 

CT'^^^^X, l^j., Srmn) = A^{x, /.^) log (^^^ + ^^^^-^(x, /X^) , (7) 

where 

Aa{x,fiF) = Ap^a (X, ftp) , (8) 

P 

= ^i?-_!:r(a:,/i^), (9) 

p 

and A'^c denotes the number of colors. The sum runs over p = q,q, g. The functions A 
and B can be expressed as convolution integrals over the parton distribution functions and 
their explicit forms are shown in Appendix^ Although Aa is scheme independent, Ba does 
depend on the mass factorization scheme, and so does the crossing function. 

After introducing the crossing function, we can write the NLO differential cross section 
in the PSS method with one cutoff scale as 

d'^mS, = J2j ^a;ic?X2|/fna:i,/XF)/f^(x2,/iF)c?af,^''(xi,X2,/i/j) (10) 

a,b 



Here da^f^^^ consists of the finite effective all-partons-in-the-final-state matrix elements, in 
which partons a and b have simply been crossed to the initial state, i.e. in which their 
momenta —pa and —pi, have been replaced by Pa and pi,, as shown in the Fig. |21 The 
difference between cia^^*^ and da^^^^'^ has been absorbed into the finite, universal crossing 
function C^{x,^f). Defining an "effective" NLO parton distribution function J^^(x) as 

^ai^) = /f + + 0(«^), (11) 

we can rewrite Eq. (jlUj) in a simple form as 

d^mS, = E / rfa;idx2^fn^i)-^f^(^2)rfa,^^^(xi,X2). (12) 

a,b 

C. 75 problem 



Because of the presence of the axial- vector current, a prescription to handle the 75 ma- 
trices in d{= 4 — 2e) dimensions has to be chosen. In this paper, we show the results 

in 



of our calculations using both the dimensional regularization (DREG) scheme ('t Hooft- 
Veltman scheme |50|) and the dimensional reduction (DRED) scheme (four dimensional 
helicity scheme j^) to regulate the ultraviolet and infrared divergencies presented in the 
NLO calculations. We note that the results of form factors and the crossing function should 
be done consistently in a given scheme. Except for the top quark mass renormalization, 
we work in the MS scheme throughout the paper to perform the needed renormalization 
and factorization procedures in order to calculate any ultraviolet and infrared finite physical 
observable. To renormalize the top quark mass, we use the on-shell subtraction scheme. 



III. LEADING ORDER RESULTS OF SINGLE TOP QUARK PRODUCTION 
AND DECAY PROCESS 



In this section we present the leading order results of the single top quark processes. 
Using the density matrix method in the NWA, cf. Eq. Q, we factorize the s-channel and 
t-channel single top quark processes (cf. Fig. H} into the top quark production and decay, 



separately. To compute the amplitudes we use the spinor helicity methods (72 



m 



75 



with the conventions as in Ref. [71|, and for completeness, we briefly review the notation in 
Appendix |X1 We note that in some of Refs. [l^, ll^ Q, 13, ll^ the phase conventions do 
not correspond to the helicity convention utilized in this paper. Below, we give the explicit 
Born level helicity amphtudes of the single top quark production and decay, respectively. 



A. Helicity matrix elements of single top quark production 

The helicity amplitudes for the s-channel single top quark production can be written as 
following: 



Mr\>^t = +) = 2 (i+\d+j (^u+\b-) ^i, (13) 
j^prod^^^ = _) = 2 (^i- \d+^ (^u+\b-^ uX, (14) 

where we have suppressed, for simplicity, the common factor \/'lEu\/lEg\/ 2Ei„ the coupling 
constants (^~^^ propagator — with s = (pu + pj)^. Here, g is the SU{2) 



coupling constant, mw denotes the mass of VT-boson, and cu^. = y-Et ± |pt |, where Et and pt 
are the energy and momentum of the top quark, respectively. The meaning of the bra (<|) 



11 



and ket (|>) in the above hehcity amphtudes is summarized in Appendix^ We note that 
u, d, t and b within the bra and ket denote the normahzed three-momentum of the particle, 
cf. Eq. ()A6jl . We did not write explicitly the helicity states of the other massless quarks 
because only one set of the helicity states give a nonvanishing matrix element. For example, 
in this case the incoming u-quark is left-handed, d is right-handed and b is right-handed. To 
calculate the squared matrix element, one also needs to include the proper spin and color 
factors which are not explicitly shown in this paper. 

For the t-channel single top quark production process, the helicity amplitudes are given 

by 

Mr'i^t = +) = 2(n+|6-) (15) 
Mr'i^t = -) = 2(u+\b-) (i-\d+)ujl, (16) 



where we have also suppressed the common factor \/ 2Eu\/ IE 2E^^ the coupling constants 
^-^^ , and the propagator ^ with t = {pu — Pdf- 

B. Helicity matrix elements of top quark decay 

For the top quark decay process, the helicity amplitude are given by 

M'''\\t = +) = -2(b'-\0+'^{e+\i+)ujt, (17) 
M'^-(Ai = -) = -2(b'-\iy+^{e+\i-)ujX, (18) 

where we only consider the leptonic decay mode of W boson and suppress, for simplicity, 
the common factor i/ZE^A/ZE^y^ZEb^, the coupling constants ^"^^ ' propagator 

5-- , where pw and Tw are the 4- momenta and the total decay width of 

ly-boson, respectively. All through out this paper we use b' to denote the bottom quark 
from top decay. 



IV. NLO MATRIX ELEMENTS OF SINGLE TOP QUARK PRODUCTION AND 
DECAY PROCESSES 



Beyond the leading order, an additional gluon can be radiated from the quark lines or 
appear as the initial parton in the single top quark process . Since the single top quark can 

19 



only be produced through the electroweak interaction in the SM, we can further separate 
the single top quark processes into smaller gauge invariant sets, even at the NLO. Taking 
advantage of this property, in the first part of this section we separate the s-channel and 
t-channel single top quark processes into smaller gauge invariant sets of diagrams to organize 
our calculations. As we pointed out in Sec. IIIBl NLO QCD corrections in the PSS method 
can be separated into two parts: (I) the resolved real emission corrections and (II) the virtual 
correction plus the unresolved real (soft+coUinear) emission corrections, denoted by "SCV" . 
After integrating out the virtual gluon and the unresolved partons, the SCV corrections can 
be written as form factors multiplying the Born level vertex. The form factors either modify 
the Born level coupling or give rise to new Lorentz structure of W coupling to fermions. 
In the second part of this section, we will write down the most general form factors of 
the single top quark processes and show their contribution to the helicity amplitudes for 
both s-channel and t-channel processes explicitly. It is worthwhile to mention that the form 
factor formalism presented here can be easily extended to study new physics models whose 
effects also show up as form factors. The derivation of the form factors for single top quark 
production and decay processes as predicted by the SM can be found in the second part of 
this section. The resolved corrections are also calculated using helicity amplitude method 
and the results are shown in the third part of this section. 



A. Categorizing the single top quark processes 

Here, we separate the NLO s-channel and t-channel single top quark processes into smaller 
gauge invariant sets of diagrams to organize our calculations. The NLO s-channel diagrams 
consist of all the virtual correction diagrams as well as the Feynman diagrams of the following 
real correction processes: 

qq' ^W*g ^bgt{^bW+), (19) 
qg W*q' bq't{^ bW+), (20) 
gq' W*q bqt{^ bW+), (21) 
qq -^W* bgt{^ bW^), (22) 
qq ^W* ^ bt{^ bW+g), (23) 



1.3 



where the gluon is connected only to (g, q') hnes in ()19 |) -(PT |) . and the gluon connected only 
to (t, fo) lines in (j^^ and the gluon connected only to (t, 6) lines in ((221) • We note that 
diagrams ()2()|) and ()21|) do not include those in which the gluon line is connected to the 
final state h and t line, for those are part of the NLO corrections to t-channel process as 
shown in Eqs. fl26p and (j2Zj)- To facilitate the presentation of our calculation, we separate 
the s-channel higher order QCD corrections (including both virtual and real corrections) 
into the following three categories: 

• corrections to the initial state of the s-channel single top quark production (INIT), in 
which the gluon is only connected to the initial state light quark (g, q') line, 

• corrections to the final state of the s-channel single top quark production (FINAL), 
in which the gluon is only connected to the final state heavy quark (t, h) line of the 
single top quark production, 

• corrections to the decay of the top quark (SDEC), in which the gluon is connected to 
the heavy quark (t, h) line of the top quark decay. 

The three types of corrections are illustrated in the upper part of Fig. |21 in which the blobs 
represent the higher order QCD corrections. The explicit real emission diagrams for the 
s-channel process can be found in Fig. |3] 

The NLO t-channel real correction processes for the top quark production and decay are 



hq — 




bW+), 


(24) 


bq' - 


Qgt{- 


bW+), 


(25) 


Qg - 


^ q'bt{- 


bW^), 


(26) 


Q'g - 


qbt{- 


bW^), 


(27) 


bg - 


^ Q(l't{- 


bW^), 


(28) 


bq — 


. q't{-. 


bW+g), 


(29) 


bq - 


-^qt{-^ 


bW^g). 


(30) 



Here the gluon is connected to both (g, q') lines and (t, b) lines in (|24[ . but only to (t, b) 
lines in ITTj) . In (|^. we restrict the gluon to be connected only to (g, g') lines. When 
the gluon in ()28|) is connected to (t, b) lines, it corresponds to the process bg tW with 
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INIT FNAL 




LIGHT 



HEAVY 



TDEC 



Figure 3: The way we organize our calculations at the NLO. The blobs in the diagrams denote the 
higher order QCD corrections, including both virtual and real emission contributions. 

W — > gg', therefore it is not included here. As done in the s-channel case, we separate the 
t-channel NLO QCD corrections (including both virtual and real corrections) (jSMSOI) into 
three categories. As illustrated in the lower part of Fig. |S1 they are: 

• the one in which the gluon is connected to the light quark (g, q') lines (LIGHT), 

• the one in which the gluon is connected to the heavy quark (t, fo) lines (HEAVY), 

• the one in which the gluon is radiated from the heavy quark (t, h) lines of the on-shell 
top quark decay processes (TDEC), 

The explicit real emission diagrams for the t-channel process can be found in Fig. |3| 
B. Form factor formalism for SCV corrections 

Below, we present the form factor formalism for each category defined in Sec. II V Al after 
including both the virtual and unresolved real corrections. 

1. INIT form factors 

Higher order QCD corrections to the diagrams labeled as INIT in Fig. |21 do not change 
the Lorentz structure of the W* — u — d coupling, therefore the most general form of the 
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initial state contribution can be rewritten as 

^7^Pl/l, (31) 

where II denotes the effective form factor that includes the higher order corrections. Denot- 
ing the helicity amplitude as M.init{\) with top quark helicity = ±1 and suppressing, 
for simplicity, the common factor \j2E^\j2Efi^2Ei, the coupling factors ^"^^ ^"^^ 

propagators ^ with s = (pu + Pj)^, we obtain the helicity amplitudes which include 

higher order corrections to the initial state of the s-channel single top quark production as 
following: 



Minit{+) = 2lL(i+\d+jiu + \h-jJ_, (32) 
MiNiA-) = 2h(t-\d+)(u+\h-)uo\, (33) 

where cu^ = \J Et±\'pt\, cf. Appendix 1X1 

Needless to say, when calculating the scattering amplitude of the single top quark pro- 
duction and decay process, cf. Eq. (^Q), up to the NLO, the decay matrix elements in this 
case is taken to be the Born level ones as given in Eqs. (|13|) and (|14p. 



2. FINAL form factors 

In the limit that the bottom quark mass is taken to be zero *, the most general W* — t — b 
coupling, labeled as FINAL in Fig. El is 

^ {i,{FI^*Pl + F,^*Pn) - ^^H^iF^PL + F^Pr) ] , (34) 
V2 I niw ) 

where the asterisk in the superscript of the form factors indicates taking its complex conju- 
gate. This is different from the coupling in Eq. (jHT|) because the top quark mass is kept in 
the calculation, and only the bottom quark mass is taken to be zero. Because the charged 
current interacts with massless quarks in the initial state, one can use the on-shell condition 
of the massless initial state quarks to rewrite Eq. (|34j) as 

^ {i,{FI^*Pl + P«) + UF2''*Pl + FtPR)} , (35) 

* Wc take the bottom quark mass to be zero throughout our calculation because (rrih/mt)^ can be ignored 
numerically. Strictly speaking, as In(mb) terms have been included in the definition of NLO PDF. 



Ifi 



where the mw has been absorbed into form factors F<^* and -F^*. Denoting the hehcity 
amphtude as AipiNALi^i^ ^t), "we obtain the hehcity amphtudes which include higher order 
corrections to the s-channel single top quark production as following: 



Mfinal{-,-) = 2F,^*(t-|(J+^(n+|6-)^^ + Ff ^J-I (36) 

Mfinal{+, -) = 2F/^* {i- \u-) (rf- \U) iot + Ft (d- 1 ^u-^ (^t- \b+) (37) 

Mfinal{-.+) = 2Ft (t+\d+^(^u+\b-'^ut+Fi'* (d-\'^u~'f(t+\b-^uX, (38) 

Mfinal{+,+) = 2Ft {i+\u-)(d-\b+)ujX + Ft (d-\ ^u-)(i+\b+)Lj'_. (39) 



As before, we have suppressed the common factor \/1Fu\/1Fd\/1Fi^ the coupling factors 



' 1 

, and the propagators it- with s = {pu + pg)'^- 



3. LIGHT form factors 

The effective form factor for u — W* — d, labeled as LIGHT in Fig. El takes the exact 
same form as W* — u — d in Eq. (jHH). Hence, the helicity amplitudes MLiGHTi^t) for the 
t-channel single top quark production are given as follows: 



Mlight{+) = 2LL(i+\d+j(^u+\b-ju;t, (40) 
MlightH = 2Ldi-\d+){u+\b-)j^, (41) 



where Li is the effective coupling induced by higher order corrections. Again, we have sup- 
pressed the common factor \/2E^y/2Edy/2E'b^ the coupling factors ^^^^ ^"^^ propa- 
gators - — 3—2- with t = {pu- PdY- 

4. HEAVY form factors 

The effective form factor for b — W* — t, labeled as HEAVY in Fig. |2l takes the exact 
same form as W* — t — b in Eq. (jnSj). Hence, the helicity amphtudes J^heavy{K, ^t) for 
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the t-channel single top quark production are given as follows: 



MHEAvri-,-) = 2Ht*{t-\d+}{u+\b-}iul-H^*{d-\ ^u-} {t-\b-} ujI, (42) 
Mheavy{+,-) = 



rL* 





(d- 


1 m- 


- < 


(d-\ 




-Ht 


(d- 




- < 


(d-\ 





i+\b-)uj\, (44) 



where H^'-^ denote the effective couplings induced by higher order corrections. Here, we 
have suppressed the common factor ^j2E^^j2Ed^/2Ei,, the coupling factors ^"^^ ^^"^ 

propagators 5- with t = (p„ — PdY- 

t - 

5. Top quark decay form factors 

The most general t — b — W coupling, labeled as DEC for both s-channel and t-channel 
processes, is 

^ bMPL + D^Pn) - b'^iD^Pn + D^P,)} , (46) 

where 2 denote the form factors which include higher order QCD corrections. Denoting 
the helicity amplitude as M.dec{K^ K'), we obtain the helicity amplitudes which include 
higher order corrections to single top quark decay process as following: 



MoEci-,-) = -2D^(^b'-\0+j{e+\i-)ui + D^p-\ ^\e-) (b' ut, (47) 

Mdec{+,-) = -2D^(b'-\0+^{e+\i+)ut+D^{i)-\ ^\e-) ($' uX, (48) 

Mdec{-,+) = -2D^' (b' + \e-^ iut + }/\e-) (b' + J^, (49) 

Mdec{+,+) = -2Df(b' + \e-\{u-\i+)ul + D^{u-\^\e-)(b'+\i+\uj'_, (50) 



where we ignore the common factor \/ 2Ee \/2Ey\/2Ey , the coupling factors ^^^^ 

the propagator 5— ■. — with pvy = Pe+ + Vv 

{pfy - m^) + imwTw 

The category SDEC (or TDEC) in Fig. |21 is obtained by convoluting the s-channel (or 

t-channel) Born level helicity amplitudes, cf. Eqs. (|13|) and (jl4j) (or Eqs. (|15|) and 

with the corresponding DEC amplitudes listed above. 
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Figure 4: Feynman diagrams of the real emission corrections to s-channel single top quark produc- 
tion. 

C. Helicity amplitudes of resolved contributions 

Here, we present the hehcity ampHtudes of resolved corrections for each category defined 
in Sec. WJ\ 

1. NLO corrections to IN IT 

The Feynman diagrams of the initial state real emission corrections are shown in 
Figs. IHa)-(f). At the NLO, a hard gluon can be radiated from the initial state quark 
line, or a quark can be generated from the gluon splitting. We separate the NLO INIT real 
emission corrections into three categories: 

INI - A : qq ^ W*g tig, including (a) and (b), 
INI - B : qg ^ W*q' tbq', including (c) and (d), 
INI — C : q g ^ W*q tbq, including (e) and (f), 

which are separately gauge invariant. Denoting the helicity amplitude as Aif}^j''"{\t), we 
calculate the helicity amplitudes for a given helicity state {Xt) of the top quark, which are 
listed as follows. 



IP 



The helicity amplitudes for INI A are : 



MU{+) = 2c.i <( ^^^^ / + \ ^^^^ '- ) (51) 



with p = Pu-Pg and q=Pd- Pg- 

The helicity amplitudes for INI B are : 



. b+\u-){d-\ ^+ '^_\t-) {d-\t-) {h^\ C ^+\u- 
Mf^,{+) = 2u;t{ / + \ /-^^ ^ } , (53) 



b+\u-) {d-\ ^+ fi\t+) (d-\t+) {b+-\ 4^ 
M^^A-) = <( ^ + ^ '- ) , (54) 



p2 q2 



with p = Pg-pd and q = Pg + Pu- 

The helicity amplitudes for INI C are : 



g2 p2 



M%A-) = ^^\{ ^ + ^ ^ ) , (56) 

with p = Pg-pu and q = Pg + Pd- 

Again in all the above equations, we have suppressed the common factor 

\ 2 1 



^j2Eu\/1E^-^2Ei, the coupling factors ( —= , and the propagator „ „ 
with piy = Pi + Pi- Here, is the coupling constant of the strong interaction. 



2. NLO corrections to FINAL 

The Feynman diagrams for NLO real emission corrections to the final state of s-channel 
top quark production process are shown in the Fig. Efg) and (h). Denoting the helicity 
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Figure 5: Feynman diagrams of the real emission corrections to the t-channel single top quark 
production 

amplitude as M-FiNAL^^t)^ then 

'u+\h-) f+ ^_\d- 



Mfinal{+) = 2iu[ 



2 J 



— mf 



u+\b-j {t+\f_\d+ 
— ml 



i+\d+) {u+\ C f + \b- 



(57) 



Mfinal{—) 

- 2uj\- 

with p = Pg + Pt and q = 
the coupling constants 
Pw =Pu+ Pd- 



u+\b-){t-\ f+ ^^\d+) {u+\b-) {t-\ f^\d+ 

2 ^ + 2mtut 



p"^ — ml 



p^ — m\ 



(58) 



^g^Ph- We again suppressed the common factor \j2Eu\J2E^^2E\ 



\V2 



and the W boson propagator 



Pw - + irnw^w 



with 
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3. NLO corrections to LIGHT 



The Feynman diagrams, that generate real emission contributions through couphng a 
gluon to the hght quark hues, are shown in Figs. a) to (f). To facihtate our calculations, 
we separate the NLO LIGHT real emission corrections into the following three categories: 

LIGHT — A : hq g'^ft, including (a) and (b), 
LIGHT — B : hq' ^ qgt, including (c) and (d), 
LIGHT — G : bg qq't. including (e) and (f). 
Denoting the helicity amplitude as -MLWHTi^t) ^ then the helicity amplitudes for LIGHT-A 



are: 



\t-)(b+\^_ f+\u-) /b+\u-)(d-\ f+ ^^\t- 




, , ,^\^- {b+\u-)(d-\f+(^_\t- 



p2 q2 



with p = Pu-Pg and q =Pd + Pg- 

The helicity amplitudes for LIGHT-B are: 



b+\u-){d-\ f+ ^\t-) (d-\t-) {h^\ 4_ fAu- 
Mf,aHT{+) = 2ut { ^ ^1-^ / + \ /^^-^ L 




MfjoHTi-) = { ^ — + ^ — ^^^^ m 

with p = p^- Pg and q = Pa + Pg- 

The helicity amplitudes for LIGHT-G are: 



d-\t-) {b+\ \ d- \ ^+ ^^\t- 



d-\t+) (b^\ ^_ ^+\u-) (b+\u-)(d-\ ^+ ^-\t- 
MZchA-) = 2^^ { / + \ M_ f_ ^(p4) 

with p = Pg -Pu and q = Pg - Pd- 



Again in all above equations, we suppressed the common factor v2^v2^v2^! the 
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coupling coustants 9. X ) , and the W boson propagator ^ wth 

V V 2 / Pw- f^w + ^T^W^W 



Pt - Pb- 

9.9. 



4. NLO corrections to HEAVY 

The Feynman diagrams, that generate real emission contributions through couphng a 
gluon to the heavy quark hues, are shown in Figs. Efg) to (n). We separate the NLO 
HEAVY real emission corrections into the following four categories: 

HEAVY — A : hq ^ Q'gt, including (g) and (h), 

HEAVY — B : bq ^ qgt, including (i) and (j), 

HEAVY -C : qg^ q'bt, including (k) and (1), 

HEAVY — D : (fg —* qbt, including (m) and (n). 

Denoting the helicity amplitude as J^HEAvvi^t)^ then the helicity amplitudes for HEAVY-A 
are: 



u+\b-) {t+\ f+ 4-\d+) (u+\b-) (t+\f4d+ 

-MWy(+) = 2a;i^ ^ ^ L + 2mtUj\-^ '-^ ^ 

q^ — mf q^ — mf 

i+\d+) (u+ \ ^_ f+\b- 



+ ^, (65) 

u+\b-) 4-\dA (u+\b-) f^d- 

MjjEAVY^-) = '^^i- —2 2 ^- + 2mtu;'_^ ^-^ 

^ q^ — mf q^ — mf 



t-|d+) (n+l ^„ f+\b- 



+ 2a;^^ '-^—^ (66) 

with p = Pb-Pg and q = Pg + Pt- 

The helicity amplitudes for HEAVY-B are: 

u+\b-) f+ 4-\d+) (u+\b-) li+\f4d 

M''heavy{+) = 2^-^ —2 2 ^- + 2mtuj\^ ^ 

q^-mf q^-mf 



t+\d+){u+ \ ^„ 

+ 2J^^ '-^—^ ^, (67) 



p2 



u+\b-) f+ 4-\d+) {u+\b-){i-\ f_\d+ 



M'^EAVvi-) = ^^i- —2 2 - + 2mtUjt 2- 



i-\d+}{u+ \ ^„ f+\b- 



+ ^ui ' • ' ' ^ (68) 



The helicity amplitudes for HEAVY-C are: 



MHEAVYi+) = 2^-^ ■ —2 2 ^-2mt^^ J- 



t+\d+)(u+ \ ^„ ^+\b- 
2ujt^ L, (69) 



p2 



u+\b-) ^+ ii_\d+) (u + \b-) \ ^_\d 



M%eavyH = ^^i- —2 2 ^--2mtujt 2 

q"^ — mf q^ — mf 



i-\d+} {u+ \ i)_ ^+\b- 



where p = Pg — Pi and q = Pg — Pt- 

The hehcity amphtudes for HEAVY-D are: 



u+\b-){t+\ ^+ {u+\b-) {t+\^4d 



q iiif q in^ 



p2 



where p = Pg — Pi and q = Pg — Pt- 



' ' ' 2 — (70) 
p^ 



t+\d+){u+\ 

2J_^ (71) 

u+\b-)li-\ ^+ 4^\d+) (u+\b-)li- \ 

-Mgi.Ayy(-) = 20.^^ ^ -2 L.2m,jJ 

q'^ — mf q^ — mf 



t-\d+}{u+ \ i)_ 

2..^^ (72) 



Again in all above equations, we suppressed the common factor \/2Eu\/2Ed\/2E\y, the 

coupling constants gg ( ^= ) , and the W boson propagator 5 — with pw = 

V V 2 y Pw- mfy + imwTw 

Pu - Pd- 
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(a) (b) 

Figure 6: Feynman diagrams of the real emission corrections to top quark decay processes 
5. NLO corrections to top quark decay 

The Feynman diagrams for NLO real emission corrections to the top quark decay process 
are shown in Fig. IHl Denoting the hehcity amphtude as AioEci^t), then 



b-\u+) (e+l 

Mdec{+) = 2u;i^ ^— 5 '- + 2mtu'' 





|z>+^ (e+ 




\t+) 




i+){b-\ 


ml 


-) 




\^+) (e+ 


1 r+ 


l*-> 




p2 _ 
i-){b-\ 







— ml 



+ 2ujI ^—^ ^, (73) 

(b-\v+^ (e+l 



Mdec{-) = 2cj+^ 2 + 2mtuj 



p^ — ml 



+ 2ujX ^-^ : (74) 



with p = Pt—Pg and q = Pb + Pg- We again suppressed the common factor y/2E^y/2E^y/2Eb 



± 



2 



the couphng constants gs ( ^= 1 , and the W boson propagator — ^ — with 

PW = Pe+ +Py 



V. NLO SCV FORM FACTORS OF THE SINGLE TOP QUARK PRODUCTION 
AND DECAY PROCESSES 

In this section the analytical results of the effective form factors are given in details 
together with the corresponding phase space boundary conditions which slice the phase 
space of real emission corrections into unresolved and resolved regions. Provided with such 
phase space boundary conditions, one can use the helicity amplitudes given in the previous 
section to perform numerical calculations. Since the unresolved regions of massless partons 



differ from the ones of massive partons, we separately consider the massless and massive 
partons and present the detailed derivations of the SCV form factors in Sec. IV Al and Sec. lV B( 
respectively. For comparison, we present our results in both the DREG and DRED schemes. 
We note that the form factors and the crossing functions should be applied consistently in 
any given scheme. 



A. NLO corrections to INIT 



Let us first examine the initial state corrections to the s-channel single top quark process, 
cf. Fig. 01 After calculating the effective matrix element with all the partons in the final 
state, we cross the relevant partons into the initial state to obtain the needed matrix element. 
In dimension d = 4 — 2e, the NLO matrix element for the vertex q — q' — W* can be written 



as 



(75) 



where u{q){v{q')) is the wave function of q{q'), Pl = {I — 75)/2. The calculation of the 
virtual corrections for the vertex q — q' — W* is straightforward and after renormalization it 
yields 



n 



qq'^W{V\vt) 



^ I 2 2 s 3 
—CpC^ { — ^ + -ln— 2 



47r2 



3 In A - In^ — + m'^wivnt) 

"^"^ ^2 m2 + -'scheme 



(76) 



where s = 2pg ■ pqi, Cp = 4/3, 



Ann' 



r(l + e), and the scheme dependent term 



qg'^iy(virt) . 
scheme 



^gg'^Vl/(virt) 
Scheme 



(77) 



-8 in DREG scheme, 
—7 in DRED scheme. 
We have neglected all the possible imaginary parts in the above result and also in what 
follows, because they do not contribute to cross sections up to the NLO. Note that the tree 
level amplitude corresponds to setting f^^^^ = 1 in Eq. (f7^. 



t It yields II = f! 



qq'^W 



in Eq. 
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Figure 7: Illustration for crossing the initial state partons into the final state in the process qq' 
W*g. 



In the phase space slicing method, the soft and collinear singularities in the virtual 
corrections, the poles of e in Eq. ()76|). should be canceled by the unresolved real emission 
corrections from processes shown in Eqs. (fT^ - (PT|) . Below, we will partition the phase space 
of the real emission corrections to calculate the unresolved contribution. 

As an example, let us examine the qq — W*g process. After crossing all the initial state 
partons of the process qq' — > W*g into the final state, the particles' momenta are assigned 
as in Fig. [7| which implies the crossed process W* — qq'g. Let us consider the whole phase 
space of the process W* qq' g as the identity and partition it into three regions as shown 
in Fig. El 



q'g\ '^■^rnin) 



1 = Q{\Sqg\ + \S 

= J-1 + ^2 + ^3, 



e(2s. 



Sqg\ 



Sq'g\) - 0(1 Vsl 

-0(1 V.I 



\Sq'g\) 

'^^^min)^(^^^min I '''99 1) 



(78) 



where 



-^1 

^3 



0(1 "^99 I ~^ I'^'j'sl 

0(|'5gg| 2Smjn)0(S}7-iM 



0(1 "^99 1 '^^min)^i^^min \^q'g\) ~^ 0(l'^f?'9l 2s,^jj^)0(s 



I '^'?9 I ) 5 



SqgD- 



(79) 
(80) 
(81) 



Here is the Heaviside step function and Sij = 2pi-pj, where Pi is the four-momentum of the 
particle i. In the phase space region constrained by function JF^ (resolved), there is no soft 
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Figure 8: The Sqg — Sqig plane for quark pair annihilation to virtual VF-boson showing the delineation 
into soft (J^2)) coUinear {J-2) and resolved region ( ^1). 

and collinear divergencies, therefore it can be calculated in four dimensions numerically. The 
soft region is defined by the function T21 which have both the soft and collinear divergencies. 
The collinear regions is defined by the function JF3 as shown in Fig. |H1 which only have the 
collinear singularities but no soft singularities. In function JF3, the first term denotes the 
collinear region oi g || q' and the second term represents the collinear region ol g \\ q. 

Under the soft approximation, i.e. in the soft region {^2)1 the squared matrix element 
can be written as a factor multiplying the squared Born matrix element: 



e(2s™„ - \sqg\ - Ivsl) \M{W' 



qq 



(82) 



where we have defined the eikonal factor f^*^ii'9 



W*-*qq'g 



gsCpH 



2e 



4(2pg ■ Pg' 



(83) 



{2pg ■ Pg){2pg> ' P g) 

It is very simple to analytically integrate the eikonal factors in d dimensions over 



the soft gluon momentum {tO] and get the soft factor 



'soft 



9t 



167r2F(l -e) \ Sr. 

r2 



2\ 2 



TT 



+4 In^ 2 + In^ 

3 



4 In 2 _ 2 / Smi 
e e \ s 



+ 41n2 In 



(84) 
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In addition to being singular in the soft gluon region, the matrix elements are also singular 
in the coUinear region (JF3) where the matrix elements exhibit an overall factorization. In 
the limit g \\ q, we define 

Pg ^^Ph, {1- OPh, (85) 

with Ph = Pg + Pq- In this limit, 

6(1 v,l - 2s™„)e(s™„ - \M{W* ^ qq'g)\' ^ \M{W* qq')? , (86) 

where the coUinear factor c^^^^ is defined as |3| 

= gyc^^ — —■ (87) 

^Pg ■ Pg 

The function P^a^? jg related to the Altarelli-Parisi splitting function, which depends on the 
regularization scheme. In this paper, we adopt two schemes: the conventional dimensional 
regularization (DREG) scheme and the dimensional reduction (DRED) scheme. We have 

in DREG scheme, 

(88) 

in DRED scheme. 

After integrating over the coUinear phase space [20] for the case g \\ q', we obtain the coUinear 
factor 

jW*^qq'g ^ 9s Cf /fvr^Y 
-1 167r2r(l-e) \s^,J 

2 Smi, ri \ 3 2,TT ^, 2 f '^^min\ tVF*— >g(j''{col) 




^ In J + ^ - ^ - 2 In^ J + 1 , (89) 

where the scheme dependent factor 1^^^^ 

- ,/ n 7 in DREG scheme, 

/sTwr^™'^= (90) 

I 6 in DRED scheme. 

Summing over the soft and coUinear factors and crossing the needed partons into the 
initial state, we get the contributions to M.'^^'^^* from the unresolved real (soft+coUinear) 
corrections from processes (fT^ -()21 |) as following: 

ff-^(-^^) = ^CfC, f 4 - - 1^ a + - - — + 2 In^ 2 
Air \e'^ e mf e 3 

,2 



3 In - 2 In^ + In^ + 4^:1'-'' ) , (91) 



9Q 



where the scheme dependent term I slheme'^'^^^^^ 



rgg'^W(real) 
Scheme 



7 in DREG scheme, 



(92) 



6 in DRED scheme. 



It is clear that the divergencies of 



qq'~*W{vat) 



and fl'' cancel with each other and the 



sum is finite and Smin dependent. The remaining unresolved real corrections for qq' —>■ W* 
are included through the process independent, but Smin and factorization scheme dependent 
universal crossing functions. 

The corrections from the resolved regions of processes (jl9p -(l2i p can be obtained by multi- 
plying the following phase space slicing functions to the corresponding phase space elements 
and matrix element squares in the cross section calculations: 



1 ©(Smiri Ngg'l) 
1 0('5miri I '^9'? I) 



for qq' — > W*g tbg, 



for qg — * W*q' — > tbq , 



for gq' —* W*q — > tbq. 



(93) 
(94) 
(95) 



The functions ensure the amplitude squares to be finite in four dimensions. Therefore, 
they can be calculated numerically. 

B. NLO corrections to FINAL 



Now we examine the final state corrections to W* — t — b. The NLO matrix element for 



the W* — t — b vertex can be written as ^ 



M 



PLv{b). 



(96) 



V2 ' ' ' ' ■ mt 
The above formula is valid only when W boson is on-shell or off-shell but coupled to massless 
quarks because we have neglected the term proportional to {pt + Ps)/^- At the tree level. 



t It yields F^* = /^'^'^ F^* = 2f^'^*'^/mt and * = F^^* = in Eq. 
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Figure 9: The stg — sig plane for quark pair annihilation to virtual T^-boson showing the delineation 
into soft (J^2)) coUinear (^2) and resolved region (^1). 

fW^tb ^ I g^^^ jw^tb ^ Q ^j^g j^Lo, the virtual corrections to f^^^~^ and f^^^^ are, 



respectively, 



fl 



W^tb{virt) 



An 
+3 In 



m 



^2+ ^Scheme 



W^tb{vivt) 



a. 
Air 



ric^C',<i^ln^,, 
s mt ' 



(97) 
(98) 



W^tb{virt) 



where s = {pt + Pi) \ si = 2pt ■ pi = s — rnf and the scheme dependent term /scheme 

—6 in DREG scheme, 



IS 



Scheme 



(99) 



-5 in DRED scheme. 



In the presence of massive top quark, the structure of collinear and soft poles of FINAL 
matrix elements is completely different from the massless case (INIT). The top quark mass 
serves as a regularizer for collinear singularities. Thus, the matrix element contain fewer 
singular structures. However, the presence of the top quark mass leads to more complicated 
phase space integrals. Again, let us consider the whole phase space of the process W* —>■ tbg 



as the identity and partition it into three regions as shown in Fig. El 

1 = ^{Stg + Slg — 2Smin) + ^i^Smin ~ Stg ~^ S^g) 



where 



'bgJ 



0('5tg O (s^jn s 



bg) 



^2 



e(2s. 



Stg 



■•bg), 



•^3 ©('''tg 2Sm,jn ) (Sffiin, Sf)g 



(100) 

(101) 
(102) 
(103) 



Here again, we divide the phase space of process W* tbg into three parts: the resolved 
region (^i), the soft region (JF2) and the colhnear region (J^a). Moreover, the phase space 
boundary conditions are much simpler than the case of massless partons, cf. Eqs. (|79| ) -(|81| ) . 

Under the soft approximation, in the soft region (^^^2), the squared matrix element can 
be written as a factor multiplying the squared Born matrix element: 



1 2 Pg^O^ fw'^tbg 



e(2s™„ - Stg - Slg) \Miw* ^ %)p /; 

where we have defined the eikonal factor as j3]: 

fW*^tbg _ „ ^ ,,2e 



M{W* ^tb)\ 



(104) 



4(2pt ■ pi) 



4m? 



X / NO ■ (105) 

{2pfPg}{2pi-pg) {2pfPgY_ 

Integrating the eikonal factors f^*^*^^ in d dimensions over the soft gluon momentum [3| ■, 
we get the soft factor 

^2 C, 



jW*^tbg 
-'soft 



_9l 

167r2r(l - e 

1 1 

e 



In ( 1 + ^ 

mi 



21n22 - 21n2 



Si + mf 
f 2 In 2 - 1 



21n2 + ii±^ 



In 1 



Sl 



mt 



Sl 



Sl + m^ 



(106) 



In the collinear region (J^s), where g \\ b, the matrix elements exhibit an overall factor- 
ization as 



e{stg - 2Sm.in)Q{s^in " %) |A^(W^* tbg) 



1 2 g\\h ig^i 



{Cj\M{W* ^tb)\\ (107) 
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where the colhnear factor c^^^^ is defined as j?^: 



:io8) 



2P5 ■ Pg {2pt ■ 

in which P^5~*^(^) is same as Eq. (|88p . Integrating over the colhnear phase space we 
get the colhnear factor 



jW*^tbg 
-'col 



Jl 

IQ-K^Vil -e) 



X 



In^ 



2s • 



m 



t , ;-Vl/*-»t53(col) 



Si 



Scheme 



(109) 



where the scheme dependent factor -^^sThemef^*'™''* 

'7 



j^W*^tbg{col) 
Scheme 



in DREG scheme, 



3 in DRED scheme. 



(110) 



Summing over the soft and colhnear factor, we get the contributions to M^*"*^^ 
the unresolved real (soft+collinear) corrections as following: 



from 



W*^tbg{rcal) 



47r 



— C^a <^ ^ + In ^ 



27r2 



2e 



mt 



.2Li2(^)-Mn^ + 21n21n 



3 

h 



+ ln^2-21n2 

,2 



mr 



mt 



(2 



2m^ , , s 
—^) In- 



mr 



ln2 j;, - In' 



2 -51 



mr 



+ 21n^ln^-^ 

mf mf Sl 



^W*^thg{rcai) 
Scheme 



;iii) 



where the scheme dependent term I 



W*^thg(rcai) 
Scheme 



is 



r VF*^%(real) 
Scheme 



- in DREG scheme, 
2 

3 in DRED scheme. 



(112) 



The correction from the resolved regions of process ()22p can be obtained by multiplying 
the following phase space slicing functions to the corresponding phase space elements and 
matrix element squares in the cross section calculations: 



0(st(; ~l~ Sjj„ 2sj7jjfi) 0(s(g 2smjn) O (Sffjjj^ s 



bg) 



;ii3) 



C. NLO corrections to LIGHT 



The NLO matrix element for the q — W* — q' vertex can be written as ^ 

Mr^*"' = ^uiq'M-'^^^'^.P.Hq), (114) 



where u{q){u{q')) is the wave function of q{q'). The virtual correction to ff^^ is 



rq^W*q'{yirt) J 2 2 , , — t , 3 

^2 



An 1 e mi e 



-i 

3 ' ^m'}' 



(115) 



where i = —2pq ■ pg/ and the scheme dependent term I^di^e '^^^'^ 

. —8 in DREG scheme, 

I —7 in DRED scheme. 

The tree level amplitude can be obtained by setting f1^^ = 1. 

We now consider the unresolved real correction to ff^^ . There are two processes that 
contribute to g — W* — q' vertex: 

• bq tq'g{I), in which the gluon only connects with the light quark {q,q') line, cf. 
Eqs. (121 and (1231), 

• bg ^ tqq', cf. Eq. 

The soft and coUinear divergent regions of bq — > tq'g (/) can be constrained by the function 

0(2Sj7iiri I '^99 1 I'^'/'sl) 

( I "^gg I 2Sj^j,^) (Smin I'^g'gl) 
+ Qi\Sq'g\ - '2Smin)Q{Smin " \Sgg\) ■ (117) 

In the above function, the first term constrains Pg to be soft and the second and third terms 
restrict Pg to be collinear with Pgi and Pq, respectively. The process bg tqq' has only 
coUinear divergent phase space region which is projected by 



0("Smm N'jsl) ©("^mm I'^g'gl) 



;ii8) 



It yields II ^ /f-*^*"?' in Eqs. ^-g^. 



in which the two terms require Pg to be colhnear with pg and pq/, respectively. After per- 
forming all the above constrained phase space integrations analytically, one can get the 
contribution to fi^^ ^ from the unresolved real emission corrections as: 

,-.W','ire.i) ^ ^^^^ ( 1 - 2 3 _ 4vr^ ^ 2 In^ 2 

4:71 ye^ e mf e 3 

+31nf^ - 21n^(-=^) + ln^(-^) + li:Z^'-''] , (119) 
where the scheme dependent term I^2!sme 



u ,s 7 in DREG scheme, 

J g— >W g (real) I ' 

Scheme i 

6 in DRED scheme. 



(120) 



It is clear that the divergencies of fl~^^ and fl^^ '''^^'^^'^ cancel with each other 

and the sum is finite and Smin dependent. The remaining unresolved real corrections for 
q — * W*q' are included through the process independent, but s^m and factorization scheme 
dependent universal crossing functions. 

The resolved phase spaces without divergent regions are obtained by multiplying the 
following function to the phase space: 

~^ I'^Q'sI 2Smin.) ^(I'^'/'jl 2Sjnin) © ("^min I'^g'sl) 

for bq^tq'g{I), (121) 
for bg tqq'. (122) 



O ( I '^(jr'g I 2Sniin) O ('^mm I "^99 1 



1 0(Smjn ["^q'sD 0('5mm I "^93 1) 



D. NLO corrections to HEAVY 



The NLO matrix element for the h — W* — t vertex can be written as ^ 



V 2 ITT-t 



The above formula is valid only when W boson is on-shell or off-shell but coupled to massless 
quarks because we have neglected the term proportional to {pt—Phj^i- At the LO, f^^*^^ = 1 



^ It yields H[* = /j*^*^*, H^* = 2/|^*^7mt and iJf* = H^* = in Eqs. 



and /a^*^* = 0. At the NLO, the virtual correction to f^w'^t /|^*^* are 

47r I f 2e e mf ti 

iQi "^t 1 1 2/~^i\ I r''W'*^*(virt)l no/1^ 

+31n(^)--ln(^)-ln(-^) + /3,,_^ /' ^^24) 

where ii = i — mf = —2ph ■ pt, and the scheme dependent term scheme* 

. —6 in DREG scheme, 

-'Scheme " \ U^Oj 

I —5 in DRED scheme. 

We now consider the unresolved real correction to f^^*-"^. The unresolved real correction 
to f^^*^^ comes from the soft and coUinear regions of the following three processes: 

• bq —>■ tq'g{II), in which the gluon only connect with the heavy quark {t,b) line, cf. 
Eqs. ^ and (H 



• qg tq'b, cf. Eq. 

• q'g qbt, cf. Eq. 

The soft and collinear divergent regions of bq tq'g {II) are sliced out by 



(127) 



The qg tq'b and q'g qbt processes both have the collinear divergent region restricted by 
©(■5m,m — l-^fogp- After integrating out the soft and collinear regions, we get the contribution 
to the form factor f^^*-'^ from the unresolved real emission corrections as: 

fi = ir'^'''^'\-+o In— - — + ln 2-21n2-2Li2 ^ ^ 

4:71 [e'^ 2e e mf 3 rnj - ti 

- J In ^ + 2 In 2 In(^) + (2 - ^) ln(l -\)-ln\l-^) 

- In2(^) - In^ ^ + 41n(^) In ^ + ^ + /sTre*^''"'^ ' (128) 

mf mf mj mj ti J 

where the scheme dependent term -^scheme*'''^'^'^^^ 

7 

,^ I - in DREG scheme, 
3 in DRED scheme. 



and the remaining unresolved real corrections for bW* — t are included through the process 
independent, but s^in and factorization scheme dependent universal crossing functions. 

The resolved phase spaces without divergent regions are obtained by multiplying the 
following function to the phase space: 

^d'^^'sl ~^ l^tgl '^^min) ^d'^tgl 2S}nm) ("^min \^bg 



for qg — > tq'b and 

q'g ^ qbt. (131) 



for bq tq'g {I I), (130) 



E. NLO corrections to the decay process t — > Wb' 



The NLO matrix element for the t — W — b' vertex can be written as ** 



Mlr^" = + ir'^'''^^^^^PnHt), (132) 



where Pr = (1 + 75)/2 and b' denotes the bottom quark from the top decay. At the tree 
level, fl^^''' = 1, ft^^' = 0. At the NLO, the virtual correction to fp^^' and ft^^''' 
are, respectively, 

A = -T-CpCe < — ^ - TT + ~ I'^ll ~ l^w) + 2Li2( ^ 



in(i - - - M + ii2::r''] , (133) 

/r^''^""^ = ^Cj.C, I ^ ln(l -Pw)], (134) 



47r [Pw ) 

where Pw = ^w/'^'ti scheme dependent term -^^s^ieme*'^''^*^ 

* w!,// .^^ I —6 in DREG scheme, 
4211^""^ = <! (135) 
-5 in DRED scheme. 

The unresolved real correction to f^"^^' is obtained by integrating out the soft and 
coUinear regions of t — Wb' g which are sliced by 



0(2Sm,j„ — \Stg\ — \Sb'g\) + ©d-Stgl 2Smin.) © (-Smin ~ I'^b'gl) 



(136) 



It yields D{ ^ ft"^^' , D§ = -2/^-*^''' /mt and Df = = in Eq. 

9,7 



After integrating over the sliced regions, we get the contribution to f^^^''' as 
fV'^''-'' = ^C,a { ^ + I - ^ ln(l - - ^ + In^ 2 - 2 In 2 
-2Li,(i^) - ^ In ^ + 21n21n(l - M 
' ^ ~ ln(2 - Pw) - \n\2 - (3w) - - (3w) 



w 



- M In — ^ - + ^Scheme \ > (137) 



where the scheme dependent term IsdiemJ^^^^'' 

(7 

^ ,,,,,, - in DREG scheme, 

/s:we'^^"^ = 2 (138) 
I 3 in DRED scheme. 

The resolved region of t — > Wh'g is obtained by multiplying the following function to the 
phase space: 



0(1 -5*3 1 + \Sb'g\ — '^Sm.in) ~ ©d-^tgl ~ '2Smin)^{Smin ^ l-^ft'gl) 



(139) 



We have checked the formulas of (|132|) - (|139p by comparing the result of NLO correction to 
T{t Wb') with Ref. 



VI. COMBINING THE PRODUCTION AND DECAY PROCESSES 

With those building blocks given in the above sections, the NLO QCD corrections to 
single top quark production and decay can be computed, keeping the full information on 
the spin configuration of the intermediate top quark state. The general differential hadronic 
cross section at NLO can be written as 

da(HiH2 YX) 

= ^ / dxidx2< fa^{xi, fiF)f^^{x2, fip) X [dao{ab — > F) + dai{ab Y)] 

a,b ^ I 

+asfa^{xi, fj,F)Cl;'''{x2, fiF, Smin)dao{ab -> Y) 
+asC^^{xi,^iF,Smin)fb^{x2,l^F)d(TQ{ab ^Y) + (xi ^ X2)|, (140) 

where da^ is the leading-order subprocess cross section, dai is the 0{as) "crossed" subpro- 
cess cross section, cf. Fig. |21 



We now consider the single top quark production subprocess ab t\hi with t\ Wph2 
and Wp — > Iv. (Here, hi and h2 stand for any single parton or multiple partons. A and p 
are the top quark spin and W boson polarization indices, respectively) In the frame work of 
NWA, the cross section can be written as 

da{ah ^ lvhih2) = ^\^M{ab txhi)M{tx ^ Wph2)M{Wp lu)]"^ 

X Sf^-^ ^— d$(a6 ^ thi)d^it Wh2)d<l>(W M, (141) 

zrritTt 2mw^w 

where Sp denotes the proper spin and color factors, (i$'s are the phase space elements 
{{2it)^S'^{P — "S^ Pi) TT ). At the LO, Fj in the above equation should be replaced 

by the Born level decay width r°(t hW). At the NLO, special cares should be taken to 



assign Ft, cf. Eq. fll49p . Tw is the W boson total decay width. 

The matrix element square can be calculated as follows. The sum over p (the polarization 
state of the W boson from top decay) is equivalent to the following replacement in M{t\ — >■ 
Wph2): 

e%,-^^u,l^{\-l^W (142) 

We denote the result by M{t\ Wh2). Decomposing M{ab txhi) and M{tx Wh2) by 

M{ah ^ txhi) = ux{Pt)MP'''', (143) 
M{tx ^ Wh2) = M'''%x{Pt), (144) 

where we have explicitly separated the on-shell top quark spinors from both the production 
and the decay matrix elements, then we have 

\^M{ab ^ txh)M{tx ^ Wph2)M{Wp ^ = \M'''''{^t + mt)MP"'\''. (145) 

In our calculations, M'^'^'^ and M^'^'^ are calculated numerically using helicity amplitude 
approach and can be easily obtained from the formulas presented in the sections IlIII and HVl 
Eqs. fll42j) and ()145|1 guarantee that the spin and angular correlations of the decay products 
are preserved. 
Denoting 

rf$LO = Sf^-^ ^ d<i>{ab thi)d^{t Wh2)d^{W iu), (146) 

2mjr^ 2'mwi-w 

c^$NLO _ \ ^^^^ _^ thi)d^{t Wh2)d^{W iu), (147) 

2mjrt 2mw^w 



where Tt = r°(t bW) + T]{t hW) and T\{t hW) is the 0{as) correction to the Born 
level decay width T^{t hW), the LO subprocess cross section is 

daoiab^luhh) = ^|Mo'^'=(|^i + mi)Mo^'-'^|W°, (148) 

where Mq^'^''^'^'^ stand for the LO amplitude. The NLO "crossed" subprocess cross section is 

+ j,2Re[Mtl'cviA + m,)Mr\Mt{i>t + m,)M^^'^)t] ci^^LO 

+ ^|M{^-(;^4 + mi)Mo^^"|2rf$NLO^ (149) 

where ^iscyi^ stand for the 0{as) amplitudes contributed from either 
soft+collinear+virtual or resolved real corrections for the production or decay pro- 
cesses. The first term is the real NLO correction from production. The second term is the 
soft+collinear+virtual correction from production. The last two terms are the corrections 
from the top quark decay. If no kinematical cut is applied, the last two terms cancel each 
other, which means there is no net correction to the cross section from the top quark decay. 
Because the virtual correction processes and the real correction processes have different 
phase spaces d^^^ and we calculate them separately using different Monte Carlo 

programs. 



VII. CONCLUSIONS 



Precision measurement of the single top quark events requires more accurate theoretical 
prediction. The fully NLO differential cross section for on-shell single top quark production 
has been calculated two years ago, but NLO corrections to the top quark decay process are 
not included, nor the effect of the top quark width. Since the top quark production and 
decay do not occur in isolation from each other, a theoretical study that includes both kinds 
of corrections and keeps the spin correlations between the final state particles is in order. 

In this paper we have presented a complete calculation of NLO QCD corrections to both 
s-channel and t-channel single top quark production and decay processes at hadron colliders. 
In our calculation the phase space slicing method with one cutoff scale is adopted because it 



takes advantage of the generalized crossing property of the NLO matrix elements to reduce 
the analytical calculations. After calculating the effective form factors with all the partons 
in the final state, we can easily cross the needed partons into the initial state to calculate the 
s-channel or t-channel single top quark cross sections. To respect the spin correlations be- 
tween the final state particles, all the amplitudes are calculated using the helicity amplitude 
method. The form factor approach is used for including the SCV (soft+collinear+virtual) 
corrections so that our results can also be used to study new physics effects that result 
in the similar form factors. To consider the top quark production with top quark de- 



cay consistently, the "modified narrow width approximation" , cf. Sec. Ill A| is adopted in 
our calculation. Our results are given in both the DREG ('t Hooft-Veltman 75) and 
DRED schemes to treat the 75 matrix in the scattering amplitudes which is important 
for predicting the distributions of final state particles. 

A preliminary study on the phenomenology of single top physics at Tevatron collider based 
on the theoretical framework presented in this paper was already presented in Ref. 
more detailed study on the phenomenology predicted by our calculations will be presented 
in sequential paper 

Notes added: While completing the writing of this paper, we noted that another article 
dealing with the same subject, but with different method, just appeared 
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Appendix A: HELICITY AMPLITUDES 



In this appendix we briefly summarize our method for calculating the helicity amplitudes. 
The method breaks down the algebra of four-dimensional Dirac spinors and matrices into 
equivalent two-dimensional ones. In the Weyl basis, Dirac spinors have the form 



(Al) 



where for fermions 



and anti-fermions 



(A=l) I 

= ^ ' (A3) 

with uj± = E ±\p\, where E and p are the energy and momentum of the fermion, respec- 
tively. Explicitly, in spherical coordinates, 

= {E, |p I sin0cos</), IpI sin0sin</), IpI cos^) (A4) 

The XA/2's are eigenvectors of the helicity operator 

h = p- a, (A5) 

where p = p/ \p\ and the eigenvalue A = 1 stands for "spin-up" fermion and A = — 1 for 
"spin-down" fermion. 

/ cos^/2 \ / -e^'^sin^/2 \ 

Xi/2 = \p+>= \ ^ , = \p->= , (A6) 

\^e*'^ sine/2 J ^ cose/2 J 

where we introduce the shorthand notations |j)±> for x±i/2- Furthermore, 

<p±\ = {\p±>y (A7) 

where the superscript denotes taking hermitian complex conjugation. Under the operation 
of charge conjugation, denoted as |;>+>, we have 

\p+ >= ia2 \p+ >*=-\p-> . (A8) 

4fi 



Similarly, 



l^-> = + I P+>, (A9) 
<P+\ = - <P-\, (AlO) 
<P-\ = + <P + \- (AH) 

Gamma matrices in the Weyl basis have the form 

where cxj are the Pauli 2x2 spin matrices. In the Weyl basis, ^ takes the form 
{po-a-p J J \7- y 

where 

l^ = {l,Ta). (A14) 
Appendix B: Ah AND Bh IN THE CROSSING FUNCTIONS 

There are four independent Ap^hx{x, fip) and B^'^'f^'^^x, fip) coefficient functions in the 
process independent, but Smin and factorization scheme dependent crossing functions, cf. 
Eq. ((7j). They are listed below, after suppressing the fip dependence. 

-2n/ ^ 



^™(^) = [^^ + 2\n{l-x)]f^{x) (Bl) 
+2 /" rf./f + l-z]+2 r d j9i^l^)-f 9i^) 



A 8 ^ , 4 (l + ^^)/<(x/z)-2/f(a;) 

- J dzfg [X/Z) 



Ag_,q{x) = - J dzff{x/z) (B3) 

= ^ f dzff{x/z) ^^^\;'^\ (B4) 
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+2 / dzf^{x/z) ln(l - z) [— ^ + l-z] (B6) 
+2h,Ml_,)WfhVM, ,B7) 

i Z 



4 /"^ 
+ - / (i^ ln(l — 2 



(1 + z')/zfi'{x/z) - 2/f (x) 



1 - z 

Bf!,,{x) = \ dzff{xlz) f + ~ ln(l - z) + 2(1 - z)], (BIO) 



= ^ f dzff{x/z)[ ^^^\ ln(l - ^) + 1], (Bll) 



2 



+2 / ci^/f (x/z) ln(l - ^) + l-z] (B13) 

+2 / dzhiil - zY-^—-- 1^^, (B14) 

= ^(y - ^) + ^ - (^) - ^ (B15) 

9 ix 1 - ^ 

= ^ dzff{x/z) f ln(l -z)- 2(1 - z)], (B17) 

= ^^'ci<(xA)[i±^^^ln(l-.) + l], (B18) 

where n/ is the flavor number, //f (x) is the parton distribution function of parton h inside 
hadron H . In the above, we have set Nc = 3. The subscript MS indicates the results in the 



(B9) 



MS DREG scheme while the subscript DRED indicates the results in the DRED scheme. 
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